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Considering two rigid conical inclusions embedded in a membrane subject to lateral tension, we 

on : 

study the membrane-mediated interaction between these inclusions that originates from the hat- 
shaped membrane deformations associated with the cones. At non-vanishing lateral tensions, the 
interaction is found to depend on the orientation of the cones with respect to the membrane plane. 
The interaction of inclusions of equal orientation is repulsive at all distances between them, while 
the inclusions of opposite orientation repel each other at small separations, but attract each other 
at larger ones. Both the repulsive and attractive forces become stronger with increasing lateral 
tension. This is different from what has been predicted on the basis of the same static model 
for the case of vanishing lateral tension. Without tension, the inclusions repel each other at all 
distances independently of their relative orientation. We conclude that lateral tension may induce 
the aggregation of conical membrane inclusions. 
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I. INTRODUCTION 



Biological membranes consist of a fluid lipid bilayer with embedded amphiphilic macromolecules such as integral 
proteins [jjj. Integral proteins are expected to be much less flexible than the lipid matrix. In a general sense, any 
molecule embedded in the membrane and differing in shape or elastic properties from the surrounding lipid molecules 
can be viewed as an inclusion. The phase behavior of inclusions in the plane of the membrane is determined by 
interactions between them. If the interaction energy is sufficiently large to compete with translational entropy, it can 
l— — '■ lead to lateral self-assembly of the inclusions. Attractive forces may result in a lateral aggregation of the inclusions, 

I ' while repulsion can give rise to a regular array with maximal spacing. 
£>• , Forces between membrane inclusions can be divided into two classes. The first class consists of the well-known 
direct interactions, namely electrostatic (for charged inclusions) and van der Waals forces. A second class comprises 
indirect interactions mediated by some kind of membrane deformation [^]|| • These interactions are determined by 
the shapes of the inclusions and the elastic parameters of the inclusions and the lipid bilayer. They can be static or 
, dynamic, in one case being due to equilibrium deformations and in the other to shape fluctuations of the membrane. 
Both types of indirect interactions have been theoretically studied for the case of zero lateral tension. The static 
0^ i interactions of inclusions affecting the membrane thickness || and of conical deformations affecting the membrane 
shape J 2[ have been dealt with. Dynamic interactions were treated for inclusions modifying the local bending moduli 
|,|§jf; including the case of rigid disks |,|,§. 

In the following we consider the static interaction between conical inclusions in the presence of lateral tension. Two 
sketches of a truncated cone embedded in the membrane are given in Figs. 1 and 2. The cone is assumed to be rigid 
and to impose a uniform slope on the surrounding membrane which returns asymptotically to the flat state at large 
O ■ distances. 

We find repulsive interaction at all values of the lateral tension if the two conical inclusions have equal orientation 
, with respect to the membrane plane. By increasing the lateral tension, the interaction is weakened at larger, but 
enhanced at smaller inclusion distances. In contrast, for opposite orientations of the inclusions in a membrane with 
non-vanishing lateral tension the sign of the interaction depends on the distance between the inclusions. At small 
separations the inclusions repel each other, while at large separations the interaction is attractive. With rising tension 
the attractive potential well deepens and moves towards smaller distances between the inclusions. 
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II. SHAPE AND ENERGY OF MEMBRANES WITH CONICAL INCLUSIONS 



We consider a membrane with two embedded conical inclusions. The cross-sections of the inclusions in the mid- 
plane of the membrane are circles of radius a. The centers of the two circles are separated by the distance R (see 
Fig. 3). 

In the absence of inclusions the membrane is assumed to be flat and to lie in the xy-plane of the Cartesian system 
of coordinates. We describe the membrane equilibrium shape produced by the inclusions by a function u(x, y) which 
determines the displacement of the membrane from the xy-plane in z direction (see Fig. 2). 

At the boundaries of the conical inclusions the displacement u is assumed to fulfill the conditions (cf. ^ ) 
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where the subscript i takes the values 1 or 2 for the hrst and the second inclusion, respectively. By fj and 4n we 
denote the polar coordinates related to the center of projection Ei of the respective inclusion on the xy-plane (see 
Fig. 3). The first equation in (1) describes the boundary of each inclusion as a circle of radius a whose center is at 
height hi above the xy-plane and which is tilted with respect to the z axes by an angle Pi in the x-direction. The 
second equation in (1) takes into account that due to the conical shape of the inclusion the membrane is attached 
to the circumference of the tilted circle with a constant angle ai. It is assumed in (1) that the contact angle ai and 
tilt angle Pi are small, on -C 1, Pi <C 1, so that we set tana^ = ai, tan Pi — Pi and neglect contributions of the order 
of magnitude of Pf determining the deviation of the inclusion projection from the circular shape. At large distances 
from the inclusions, r.i 3> R, the membrane remains flat, so that Vti — * for ri — > oo. 

The inclusions characterized by small on and Pi can produce only a weak deformation of the initially flat membrane, 
which means that the gradient of the function u{x,y) remains small, |Vu| 1, everywhere along the membrane. The 
membrane energy fl9| can then be written in the approximate form 



G = 



](Au) 2 + RK + -(Vit) 2 ) d 2 r 



(2) 



where K denotes the bending rigidity, 7 the lateral tension, K the Gaussian curvature and k the modulus of the 
Gaussian curvature. In our approximation, the Laplacian Au equals the sum of the principal curvatures of the 
membrane J, while |(Vu) 2 gives the increase of membrane area per unit projected area due to membrane tilt Vit. 
The integration of (2) is performed over the projected area. The membrane shape is determined by the Euler-Lagrange 
equation following from (j^) 



AAu = -Au 

K 



(3) 



We derive the interaction energy of two conical inclusions in two steps. First, we solve the shape equation (|^) 
accounting for the boundary conditions (|l]) and the asymptotic boundary condition Vu — * for — ► 00. Second, 
inserting the obtained function u(x,y) into we determine the membrane energy. Throughout this calculation we 
assume the inclusion distance R to be large compared to the radius of inclusion a and retain only the leading terms 
in a/R. 

It is important to note that the interaction energy cannot depend on the modulus of Gaussian curvature R. According 
to the theorem of Gauss-Bonnet an integral of the Gaussian curvature K over a surface is equal to the negative sum 
of the line integrals of the geodetical curvature k g over the surface boundaries (apart from a constant that depends 
only on the genus of the surface). The value of the geodetical curvature at the inclusion boundaries is completely 
determined by the radius a and contact angles ai {\k g \ — 1/a-cosoLi) and does not depend on the distance R between 
the inclusions. Therefore, the integral of the Gaussian curvature K over the membrane must be independent of the 
distance R and does not contribute to the interaction potential. 

For any given distance R between the inclusions the energy has to be minimized with respect to the heights hi 
and tilt angles Pi. This results in conditions of zero vertical force and zero torque acting on each inclusion. The two 
conditions are expressed by the equations (see Appendix A) 
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respectively. The integration in (^) is performed over the boundary of each inclusion. 

III. INTERACTION IN ABSENCE OF LATERAL TENSION 

We first consider the important limiting case of zero lateral tension, 7 = 0. The shape equation (|J) then reads 

AAu = . (5) 

To solve this equation satisfying the boundary conditions (|l|) and the conditions of equilibrium (^) we use the following 
ansatz. We consider the function u(x, y) describing the shape of the membrane in the form 

u = ui(r 1 ,(j>i) + u 2 (r 2 ,0 2 ) , (6) 

where Ti,4>i denote polar coordinates with respect to the center of inclusion i. The relationships between the polar 
coordinates related to the first and the second inclusion are (see Fig. 3) 



ri = ^/R 2 + rf - 2Rr 2 cos <j> 2 (7) 
r 2 cos0 2 -i? 

cos 0i = — ^^^^^=^^^^^= (8) 

^JR 2 + rf - 2Rr 2 cos (j) 2 

The f unct ions u\ and u 2 in (|^) are general solutions of the shape equation (|^) in polar coordinates. They are obtained 
from ( |B2| ) derived in Appendix B and have the form: 

I X \ 4. , (*) r 1 (0 j . (i) -1 , . (i) COS 0, 

Ui( r ij W = const. + Cq mri + c\ cos </>j + c 2 Inr; cos 0j + C3 



- W rrr.2.4 I r (') cos2 ^ 1 (i) COST^j (i) cosnfc 

-4 cos/^-t-c s 2 +...+c 2 „ 2 -i-c 2n+1 i^yj 
r ,- r„- 



Terms of (B2) proportional to sin nip are omitted in (J9|) because of the mirror symmetry of the system with respect 
to the xz-plane (see Fig. 3), and terms exhibiting higher than logarithmical divergence for r; — - > 00 are left out since 
they violate the boundary condition of an asymptotically flat membrane Vu — * for ri — ► 00. The only exceptions 
are rj lnr^ cos (fn and cos 0j. From the boundary condition of an asymptotically flat membrane it can be concluded 
immediately that the coefficient c 2 ^ must be equal to —c 2 2 \ The sum of the corresponding terms then diverges only 
logarithmically for n — > 00. The terms r% cos <f>i are proportional to the Cartesian coordinate x and thus describe 
rotations of the membrane as a whole pXfl . Any such rotations must be equal but opposite, i.e. c^p — — c± , to satisfy 
the boundary conditions at infinity, so that we can as well drop these terms. 

The coefficients cj 1 ' in (^) are determined from the boundary conditions (Q) and equilibrium conditions (Q) . Consider 
these conditions at the circumference of inclusion 2. To apply them we have to express the membrane shape (^) in 
the vicinity of the inclusion. The function u 2 is simply given by (^) with i = 2. To present the function u± in a 
convenient form we take (^) with i — 1 and insert (Q) and (|J) into it. In the vicinity of the second inclusion the value 
of r 2 is close to the inclusion radius r 2 sa a. Using the assumption a <C R and, consequently, f2«JJ we perform a 
Taylor expansion about the center of the inclusion projection E 2 



Ul 



= const. + c { 1] ^ln R - ^ cos <j) 2 - cos 2</> 2 ^ - c { p f i + cos 2 

• ' 1 1 ^ -iZlniZ + (1 + In R)r 2 cos 2 - ^| + -^(cos 3</» 2 - 3 cos <fe) 

LR 12R Z 



+4- ( 1 - ^(1 - cos20 2 ) + ^ - c- - + ^ cos 02 + ^ + O ^ ) (10) 



3 



The resulting expression for the membrane shape, u = u\ + Ui is 



where 



u 


_ f(2) 
— ^0 


(r 2 ) + 


A (2) (r 2 ) 


r'2«-R 








f (2) 

Jo — 


const 


+ <£ 


lnr2 — 


/! (2) = 

J 1 


(2) 
Co T9 


In r2 - 


f"2 


f (2) 


cf + 


(2) 

c\ 
»2 


- (-4 15 


f (2) 


c (2) 
— — 1 

^2 


c (2) 

_ 7 

rf 


1 L 2 '2 
12i? 2 



(11) 
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Inserting (11) into the boundary conditions (jl]) and equilibrium conditions (Q) at the inclusion 2 we obtain a series of 



(i) 

equations for the coefficients Cj . 

To account for the boundary and equilibrium conditions at the inclusion 1, we perform the same procedure as 
described above to obtain identical equations in which the index 2 is replaced by 1 and vice versa. 

The equations obtained for the coefficients Cj can be solved order by order in the small parameter a/R. The 
solutions are 



„(i) 



ctia + O 



(2) 

and equivalent results for c^ , the remaining coefficients being of third or higher order in a/R. 

We are now in a position to compute the energy of the membrane. Omitting the contribution of the integral of the 
Gaussian curvature, which is independent of the distance R between the inclusions (see above) we obtain from (S) 
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G(R)= j^J 2 d 2 



where the curvature J is given by 



J = Au = -4c 
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In the first non- vanishing order in a/R, the energy of interaction of the inclusions is 




(16) 



According to ([16]), the energy is positive and decays monotonically at all values of the contact angles a>i, ct^ and all 
distances between the inclusions R. This means that in the case of zero lateral tension the interaction between the 
rigid conical inclusions is always repulsive. The result (^) is in agreement with an earlier calculation O which in 
addition predicts a contribution proportional to R,, the modulus of Gaussian curvature. We think that there should 
be no R term (see end of Section II). 



IV. INTERACTION IN PRESENCE OF LATERAL TENSION 



We now extend the methods of the previous section to analyze the interactions of inclusions embedded in a membrane 
subject to non-vanishing lateral tension 7. The shape equation (||) can be written as 
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AAu = £ 2 Au, (17) 

where £ = y/j/n has the dimension of a reciprocal length. To find a solution of the shape equation satisfying the 
boundary conditions (Q) and equilibrium conditions (^) we use, as in the previous section, ansatz (^|) with the functions 
Ui{ti, 4>i) being general solutions of the shape equation ( |Tt[ ) in polar coordinates. These functions are taken from (B3) 
derived in Appendix B and have the form: 



Ui = const. + c$ Ko{£ri) + Cj*Vj cos^i + c$Ki(£ri) cosfa +c^ ( ^- 



-cfK 2 {in) cos20 4 + c »2^1 + . . . + c%lK n (Zn) cosnfc + 42+i £ ^ i + . . . (18) 



The coefficients of all terms of (|B3| ) proportional to sin n<f> are taken equal to zero because of the symmetry of the 
system. Also, the coefficients of terms violating the boundary condition of asymptotically vanishing gradient of the 
displacement, Vu — * for r,; — > oo, must be zero. As we do not consider rotations of the system we set 4 = 0. 



In addition, we omit in ( |18|) the logarithmic term of (B3), as the related change of the area of the membrane would 
result in an infinite energy of the lateral tension 7. 

Equation (|l8| ) transforms into (||) in the limit of vanishing lateral tension, 7^0 (i.e. £ — * 0). This can be shown 
by inserting into (|l8|) the approximative expressions of the Bessel functions K n (x) for small arguments x 



K (x)^-lnx, K x {x)^-, K n (x)*^^(-) n -^^(-) n 2 forn>2 (19) 
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The coefficients Cj in ( p.8| ) are determined by the boundary conditions (|lj) and equilibrium conditions (|4|) in a way 
similar to that described in the preceding section. For example, we present u\ in the vicinity of the inclusion 2 by 
inserting (7) and (8) and obtain after an expansion in the small parameter r 2 /R 
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We have to stress that in this case the expansions up to the second order in r 2 /R are sufficient only if £a < 1. For 
£a ^> 1, which is equivalent to the condition of strong lateral tension, the series (|20| ) converges too slowly to be 
approximated by the sum of just a few Taylor terms. This can be seen from the asymptotic expansion of the functions 
K n (x) which for large arguments are proportional to exp(— x)/s/x irrespectively of n. 

Inserting the sum (|lj) with i = 2 and (^0|) into the boundary conditions (P and the equilibrium conditions (Q) at 
inclusion 2, we obtain for a <C R and £a < 1 a set of linear equations for the coefficients cj- . Applying the same 
procedure to satisfy the boundary and equilibrium conditions at the inclusion 1 leads to an analogous set of equations. 
Solving all these equations for a <C R and £a < 1 we obtain 

4 X) = -a x a + . . . , 4 1} = - i« 2 a(£a) 2 ifi (£R) + ■■■, (21) 

4 1} = -a 2 a(Za) 2 K 2 (ZR) + ... 

and corresponding results for cjj j 4 and cf \ The coefficients c 3 ^ and Cg are given by the relations 



5 



(22) 



The interaction energy G(R) of the inclusions is obtained by integration of ^k(Au) 2 over the xy-plane. Transforming 
the area integrals into line integrals over the boundaries of the inclusions as shown in Appendix C we find the following 
dominant terms for small a/R and £a < 1 



G(R) = 2T7Ka 1 a 2 (^a) 2 K {^R) + nn(al + a^af K%(£R) 



(23) 



In the limit of vanishing tension, 7 — > (i.e. £ — > 0), this expression for the interaction energy coincides with (|l6|), as 
can be seen by expanding the Bessel functions for small arguments according to (|l9|). 

The interaction energy ( p3|) depends on the relative orientation of the conical inclusions. If the cones are oriented 
in the same direction, their contact angles a\, a 2 have the same sign. In this case the energy of interaction is always 
positive and decreases with increasing distance R. The repulsive potential is illustrated in Fig. 4 for two identical 
and equally oriented inclusions at different lateral tensions. By increasing the tension the interaction is weakened at 
large, but enhanced at small inclusion distances. 

If the inclusions are oppositely oriented, the contact angles a,\ and a 2 have different signs and the energy of 
interaction behaves non-monotonically, as illustrated in Fig. 5. The energy has a minimum at a finite separation R* 
of the inclusions. This means that the forces between the inclusions change from repulsive to attractive depending 
on the distance R. For distances shorter than R* the inclusions repel each other, while for separations R > R* the 
interaction is attractive. With rising lateral tension the separation of zero force, R*, moves towards smaller values 
and the associated potential well deepens (see Fig. 5). 



V. CONCLUSION 



To summarize, we derived an interaction energy between two conical inclusions embedded in a fluid membrane 
subject to lateral tension. For this purpose, we calculated the equilibrium shape of an almost flat membrane and 
its bending energy in the presence of inclusions as a function of their distance. In contrast to the case of vanishing 
tension, this interaction depends on the orientations of the inclusions with respect to the membrane plane. For 
oppositely oriented inclusions the interaction changes from repulsive to attractive as the separation increases, while 
equally oriented inclusions repel each other at all distances. This is very different from the case of vanishing lateral 
tension where the interaction of conical inclusions is always repulsive, independently of relative orientation. 

We did not consider in this study the contribution of thermal undulations of the membrane to the interaction 
between the inclusions dJ^,^. In the case of non- vanishing lateral tension this may be partially justified by the 
fact that the undulations are diminished by the tension. Moreover, others have found for the case of zero tension 
that the static part of the interaction exceeds the dynamic one for k(o? + a\) > 3fcT § or 1.5 kT |,| where k is 
Boltzmann's constant and T is temperature. If the tension-induced forces dominate, they should lead to interesting 
phase behavior of embedded inclusions. For example, the attractive interaction between pairs of oppositely oriented 
conical inclusions may favor the formation of clusters with a regular structure where inclusions with different signs of 
the contact angles alternate. For an estimate of the attractive interaction one may use Fig. 5. With k = 1 • 10~ 19 J 
(typical of lipid bilayers), a = 0.5 (26.8°) and £a = 0.4, the minimum of the interaction energy G(R) — V(R) ■ a\n 
is roughly —4 • 10 -21 J (~ kT at room temperature). Because of £ = \/"//k the lateral tension needed to produce 
£a = 0.4 is given by 7 = 0.16 n/a 2 . For a = 4nm and n — 1 ■ 10~ 19 J, one finds 7 = lmN/m which is below the 
known tension of lipid bilayer rupture p"l| ] . 

While our results are intuitively appealing and may be obtainable more directly, we performed a complete per- 
turbation calculation to make sure that no terms are missed. The boundary and equilibrium conditions for the 
membrane with conical inclusions are central to our calculations. They resulted in a set of linear equations for the 
coefficients of two superimposed expansions, one for either inclusion. This method is extendable to larger number of 
inclusions by using similar sets of boundary conditions. In computer-aided calculation the shape of a membrane could 
be determined with any desired precision and for any number of inclusions. 
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APPENDIX A: DERIVATION OF THE EQUILIBRIUM CONDITIONS (4) 



To derive the equilibrium conditions (4) we study a variation 

v(r, (f>, e) = u(r, (f>) + e ■ Su(r, (j>) 



(Al) 



of the equilibrium membrane displacement u{r, <fi) on a circular ring S: a < r < b, < <fi < 2n around a conical 
inclusion. To simplify the notation we leave out indices of the polar coordinates r, (f>. The variation is restricted by 
the boundary conditions (1) of the inclusion. So 



Su 
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8c + 8(3 a cos </> 
8(3 cos 4> 



(A2) 
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where e • 8c and e • 8(3 denote the changes of the hight of the inclusion center and the tilt angle, respectively. At r = b 
we set 
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d8u 
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= 0. 
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Omitting Gaussian curvature the membrane energy (2) can be written as 
G=J s (^Avf + l(Vvf) d 2 r 
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In equilibrium the energy G is minimal. So 
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By partial integrations we obtain 
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and, inserting f(v,v r ,v^,v rr ,v^,p,r) as defined in ([Ajf), are led to: 

2tt b 2n 
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The equilibrium displacement u(r, 4>) fulfills the shape equation ( 3) o f a tense membrane. So the integrand of the area 
integral in ( |A8|) is zero at e — 0. Taking into account (A2) and (A3) we conclude 



8G = 



dG 
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e=0 
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r— yfu — kAu) {8c + 8(3 a cos 0) + nrAu 8(3 cos < 
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Since SP and 5c are independent of each other we arrive at the equations 
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SG(5c) = 6c 



a— (71* — kAu) 
or 



d(f) = 



(A10) 



5G(5j3) =5(3 / cos0 



a 2 — h/u — kAu) + naAu 
or 



d<P = Q 



(All) 



which state that the vertical force and the torque, respectively, acting on the inclusion must be zero in equilibrium. 

APPENDIX B: GENERAL SOLUTIONS OF THE SHAPE EQUATIONS IN POLAR COORDINATES 

In this appendix we derive the general solution of the shape equation AAu = £ 2 Ait in polar coordinates. We 
perform the calculation in two steps. We first look for the solution J(r, (b) of an intermediate equation A J = £ 2 J and 
then solve the equation Au = J(r, </>). General solutions of the last equation are also general solutions of the shape 
equations. 

Below we consider separately the case of £ = 0, corresponding to the vanishing lateral tension, 7 = 0, and the case 
of non- vanishing £. 



1. Vanishing lateral tension, £ = 

In this case the shape equation has the form AAm = 0. A solution of the intermediate Laplace equation AJ = 
on a circular ring can be found by the method of separation of variables and reads fl2| , |l3|| : 



J(r, (j)) — ao + bo In r + (a n cos rxj) + b„ sin 



rub) r n + 



(c n cos nip + d n sin rub) r n 



(Bl) 



The general solution of the linear inhomogeneous equation Au = J(r, cj>) is the su m o f a special solution and the 
general solution of the homogeneous equation Au = 0, the latter having the form of (Bl). We obtain 



u(r, (jy) = A a r 2 + B r 2 (In r — 1) + (A± cos <j) + B\ sin <f>)r In r 

00 00 
+ (A n cos ruj) + B n sin n$) r~ n+2 + (C„ cos ruf> + D n sin ruf>) r n+2 

n=2 n=l 

00 00 

+A Q + B In r + Yj (A n cos n(j> + B n sin n<fr) r~" + (C n cos rub + D n sin 



no) r 



(B2) 



The terms with unbarred coefficients belong to the special solution, which can be directly checked by its insertion 
into Am = J. The Ao-term of (B2) corresponds to the ao-term of (Bl) etc. The terms with barred coefficients give 
the general solution of Au = in analogy to (Bl). 



2. Non- vanishing lateral tension, £ 7^ 



By applying the method of separation of variables described in [|12|,|13| also to the case of non-vanishing lateral 
tension we find the following general solution of the intermediate equation A J = £ 2 J 

00 00 
J(r, (j>) = ao-ftTo(£r) + (a n cosntfi + b n s'mruf)) K n (£r) + (c n cosn</> + d n sin 710) I n (£r) , 

n— 1 n— 1 

where /„ and K n denote modified Bessel functions. A general solution of the equation Au = J(r, <fi) again consists of 
the sum of a special solution and the general solution (Bl) of the Laplace equation Au = 0. It can be written in the 
form 
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u{r, (f>) = A K (£r) + B I (£R) + A Q + B lnr 

oo oo 

+ 2_. (An cos n( t> + B n sin ruj>) K n (£r) + (C n cos ncf) + D n sin n<p) I n (£r) 



n—1 n—1 

oo oo 



+ (A n cos n</> + B n sin n<^>) r n + (C„ cos n0 + D n sin n^>) r" , (B3) 

n—1 n—1 

taking into account that Ko(£r), K n (£r) cos n<j>, Io(£r), I n (£r) cos ncj) and the corresponding terms containing sin n<j) 
are eigenf unctions of the Laplace operator. 



APPENDIX C: REDUCING AREA INTEGRALS TO LINE INTEGRALS IN THE CALCULATION OF 

THE MEMBRANE ENERGY 



In the calculation of the energy of the tense membrane 

G = 

we encounter, due to our ansatz ([j"8|), integrals of the form 

I(f,9) = 



d 2 r 



^AfAg+lvf-Vg 



(CI) 



(C2) 



R'/E 1 UE 2 



where either / obeys A/ = j/k ■ /, which is true for the terms of ( |Tq ) containing a Bessel function, or g is a solution 
of Ag — 0, or both (see Appendix B). Ei denotes the projection of the inclusion i into the xy-plane (see Fig. 3). 
Applying a theorem of Green we may write 



J (V/ ■ Vg) d 2 r = - 

R 2 /E 1 UE 2 R 2 /E 1 UE 2 

if Tif{dg/dri) goes to zero for — * oo. Since 



2tt 



fAg d 2 r- f 



dg 
dn 



2tt 



a d<pi 

r\ —a 



■ dg 

dr 2 



a d<pi 

r 2 =a 



(C3) 



^AfAg-lfAg 



d 2 r = 



(C4) 



J? 2 /SiU£ 2 

for A/ = j/k ■ f or Ag = 0, we find 

7 



2tt 



2tv 



i I f d 9 



a ct(p2 

r 2 =a 



(C5) 
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Figure Captions: 



Figure 1: 

Conical inclusion in a bilayer membrane 



Figure 2: 

Idealization of a conical inclusion as rigid disc of height hi and tilt angle fa making a uniform contact angle on with 
the surrounding membrane. The cross section contains the axis of the cone. 



Figure 3: 

xy-plane with inclusion projections E\ and E 2 



Figure 4: 

Dimensionless interaction potential V(R) = G(R) / (a\n) of two equally oriented inclusions 
(ai = a 2 ) as a function of the dimensionless distance R/a for £a = 0.05, 0.1, 0.2 



Figure 5: 

Dimensionless interaction potential V(R) = G(R)/(o-\k) of two oppositely oriented inclusions (a\ = — 012) as a 
function of the dimensionless distance R/a. The potential well deepens with increasing £a = 0.1, 0.2, 0.3, 0.4, 0.5. 
Note that R/a — 2 means two discs in contact. Near this value the results can only be regarded as estimates. 
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